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PREFACE 



The study of the effect of resonant expansion systems on the 
propagation of long waves is presented in this thesis as a design 
problem in coastal hydraulics. Two previous studies which have been 
made on similar systems are presented in Chapter I to provide back- 
ground material on the subject. The first study presented was made 
by Lamb (1916); and it demonstrates that reflection of long waves can 
be obtained by an expansion system and that the amount of wave energy 
that is reflected is a function of the geometry of the system. The 
second study presented was made by Valembois (1953) and is based 
on a hydrodynamic impedance theory. An important point to note in 
Chapter I is that both of these studies use a scalar equation of pressure 
continuity and a vector equation for the conservation of mass to obtain 
a solution. A new solution is derived in Chapter 2 by formulating a 
boundary value problem which incorporates the same type of boundary 
condition equations noted earlier. The problem is formulated in con- 
siderable detail to demonstrate the type of problem solution technique 
which is required to solve problems arising in coastal hydraulics. The 
details given in the problem statement become important later when 
the effect of geometrical changes of the resonant system are evaluated. 
Chapter 3 discusses the inte rdependance required between Civil 



Engineering and Oceanography to effect a complete solution to a design 
problem in coastal hydraulics. Several methods for presenting wave 
data are given, and one of these methods is selected for the design 
problem being considered. Finally, a linear analysis is employed in 
Chapter 4 to evaluate the effect of the resonator on the propagation of 
long waves by combining the solution derived in Chapter 2 with the 
spectral method of wave data presentation selected in Chapter 3. 

The results of the linear analysis are presented in graphical form 
as a measure of the amount of wave energy that is reflected and 
transmitted by the resonator. 



ABSTRACT 



The response of a single resonant expansion chamber to the 
periodic pressure fluctuations of ocean waves is to act as a rigid 
vertical reflecting surface. The effective frequency band-width over 
which reflection is significant may be determined for a given harbor 
geometry by means of a graphical linear analysis incorporating the 
reflection or transmission coefficient and wave data presented by a 
Bretschneider Power Sprectrum. The most critical dimension of 
the resonator was determined to be the length of the expansion, i , 
measured transverse to the centerline of the channel; the effect of 
the resonator was maximized when this dimension was either one- 
quarter or three-quarters of the design wavelength, L,. The influ- 
ence of the width of the resonator, b , measured parallel to the 

R 

centerline of the channel was determined to be maximal when this 
dimension was one-half of the design wavelength, L f . The effective- 
ness of the resonator is also a function of the width of the main 
channel, b^, and this parameter was incorporated in the expressions 
for the coefficients. 

Because of its analogy to electrical band- stop filters, the 
possibility for improving the effect of the resonator by constructing 
several resonators in cascade appears good and offers an opportunity 
for extending the results obtained here. 



1. INTRODUC TION 



As the logistical requirements of nations grow and expand, new 
techniques for designing harbors will be required to accomodate the 
expansion of littoral ports of communications. The port areas which 
are sheltered naturally from the ocean forces and which require little 
or no additional protective structures have already been exploited. 
Designers, therefore, must address themselves to the problems of 
the less protected areas where this expansion must occur. One of 
the major concerns in an unsheltered area is the design of breakwaters. 

The purpose of breakwaters is to cause a reduction of wave 
heights in its lee (Wiegel, 1964). The conventional methods used for 
sheltering the entrance of a harbor between breakwaters may be 
divided into three general types (Figure 1. 1). All three methods 
shown require shipping vessels to maneuver in order to enter or to 
leave the harbor. This required maneuvering is both expensive and 
inconvenient to the user. In addition, the third method shown requires 
a deep shipping channel close to the shore, and may require the addi- 
tional expense of periodic dredging. The alternatives which have been 
considered to provide a direct access to the shipping channels from 
the harbor have included pneumatic and hydraulic devices as well as 
resonant expansion chambers (Wiegel, 1964). The development of the 



latter is considered here. 
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Figure l.l. Three conventional methods used to protect 
breakwater entrances. 
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Wave Reflection by an Abrupt Expansion 



Lamb (1916) investigated the partial reflection of a wave 
resulting from an abrupt change in the cross section of a channel. 

His results were based on a linear solution to the wave equation. 

(A linear solution is presented in Chapter 2). 

A right-handed cartesian coordinate system (z-axis positive 
upwards) is oriented along the centerline of a channel at the mean- 
water surface with the origin located at the discontinuity (Figure 1. 2). 
The wave amplitudes and relative velocities to the left of the discon- 
tinuity are given by 

T, = F(t- -^) + £ (t + -r- ) (l. l) 

1 c i c i 



EC = ' [F(t - — ) - f (t + — )] 

1 C 1 C 1 C 1 

and to the right of the discontinuity by 

r\ z = 0( t - 7") 

L 2 



( 1 . 2 ) 



(1.3) 



*z--t ■ ; <‘- 4 > 

2 2 

where F is an arbitrary periodic function representing the incident 
wave; f and 0 are arbitrary periodic functions representing the 
reflected and transmitted wave, respectively; and c^ and c^ are 
the wave celerities of the wave systems to the left and to the right of 
the discontinuity, respectively. The linear superposition theory is 
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Figure 1.2. Partial wave reflection by an abrupt expansion 
in a channel. 
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valid provided that the small amplitude assumption holds; i. e. , 



i 

— « 1 for i = 1, 2. 

h. 

l 



(1. 5) 



The assumption that the waves are propagating in shallow water 
requires that 




( 1 . 6 ) 



where L is the wavelength in deepwater (Kinsman, 1965). By invoking 
d 

the shallow water wave assumption given in Eq. (1. 6), the wave 
celerities may be approximated by 



c . = g h. for i = 1 , 2. 

There are two boundary conditions which must be satisfied at the 
origin (x=y=0). The conservation of mass requires that 



(1.7) 



B l h i r i= B 2 h 2“2 ; (1 ' 8) 

where B^, B^ are the breadths of the channel measured at the mean- 
water surface, and h^, h^ are the mean-water depths. The assumption 
that the fluid motion will be sensibly uniform along and parallel to the 
centerline of the channel for small distanpes (compared to a wavelength) on 
either side of the discontinuity requires that there be no sensible 
change in the mean-water surface across the discontinuity. This con- 
dition requires that the continuity of pressure across the abrupt change 
be 
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^1 = V 



Substituting Eqs. (1.2) and (1.4) into (1.8) gives 

Vl B 2 h 2 

g — [F (t) - f(t)] = g — . <f> (t) at x = 0 ; 

C 1 2 



(1. 9) 



( 1 . 10 ) 



and the substitution of Eqs. (1. 1) and (1. 3) into (1. 9) yields 

F (t) + f (t) = 0 (t) at x = 0. (1. 11) 

The ratios of the reflected wave and of the transmitted wave to the 
incident wave determine the reflection and transmission coefficients, 
respectively; i. e. , 



K = 

R F 



and K = 

T F 



These coefficients are determined by the simultaneous solution of Eqs. 
( 1 . 10) and ( 1 . 11). By means of Eq. ( 1 . 7), the reflection and trans - 
mission coefficients, respectively, are 



B 1 C 1 - B 2 C 2 
K R B,c, + Be 



( 1 . 12 ) 



and 



K 



1 1 



2 Vi 



2 2 



Vl + B 2 C 2 



(1. 13) 



The energy contained in the reflected and in the transmitted waves is 
equal to the energy of the incident wave if the sum of the squares of the 
two coefficients is equal to unity; i. e. , 



2 2 

K + K = 1. 
R T 



(1. 14) 



Equations (1. 12) and (1. 13) establish that the effectiveness of the 
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abrupt expansion for reflecting waves is dependent on the geometry of 
the channel. 

A similar procedure to the above approximation given by Lamb 
(1916) could be employed to evaluate the reflection and transmission 
resulting from an abrupt contraction of a channel and the two solutions 
superimposed in a limiting case. However, the basic assumption that 
uniform conditions exist parallel to the centerline of the channel would 
now introduce a greater error of approximation, as the dimension of 
the length of the expansion measured parallel to the centerline of the 
channel approached small multiples of the incident wavelength. These 
results obtained by Lamb (1916) indicate that substantial reflection 
occurs from an abrupt expansion in a channel and that a further investi- 
gation of an expansion system is justified. 

The Resonant Expansion Chamber 

The resonant expansion chamber (resonator) takes advantage 
of the periodic character of the surface disturbances and of the pres- 
sure fluctuations produced by a wave (Valembois, 1953). The solution 
derived by Valembois includes a wave system transverse to the center- 
line of the channel which is induced in the resonator; this wave system 
was neglected by Lamb in the preceding example. The derivation again 
assumes a linear solution to the wave equation but provides an additional 
parameter for evaluating any energy dissipation which may result. The 
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measure of energy dissipation utilizes the Neumann impedance theory 
for a hydrodynamic oscillatory system; this theory is analogous to 
the impedance theory for electrical oscillatory systems (Defant, 1961). 

The hydraulic impedance of the discontinuity (Figure 1.3) is 
in accordance with Ohm's law for the analogous electrical oscillatory 
system. The effect of the resonator is to alternately withdraw from 
and discharge into the main channel a quantity of fluid which is pro- 
portional to the surface distortion associated with the incident wave 
system. The constant of proportionality is a measure of the hydraulic 
impedance of the resonator; i. e. , 

h a 

z _ — (1.15) 

V 

A right-handed cartesian coordinate system (z-axis positive upwards) 
is oriented at the orthogonal intersection of the centerlines of the main 
channel and of the resonator at the mean-water depth. The linear 
shallow water approximation for the wave celerity may be written as 
in the preceding example as 

c = \Tgh ; (1. 16) 

where h is the mean-water depth of the channel. The relative velocity 
U(x, t) and superelevation H(x, t) in the section to the left of the dis- 
continuity at A are given by the following relationships: 



u (x,t) = F (t - — ) - f (t + — ) 

c c 



H (x, t) * -S- = F (t - — ) + f (t + — ) , 



(1. 17) 



c 



(1. 18) 
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Figure 1.3 



Action of a resonator on the propagation of a 
gravity wave, (after Valembois, 1953) 
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and to the right of A by the following: 



x 



tr (x, t) = f' (t - 



c 



(1. 19) 



H" (x, t)- = F' (t - ^-) ; 



( 1 . 20 ) 



where F is an arbitrary periodic function representing the incident 
wave; and f, F' are arbitrary periodic functions representing the 
reflected and transmitted waves, respectively. Since the medium of 
propagation is assumed to be at rest, all relative velocities will denote 
absolute velocities. The channel is also assumed to be infinite in 
extent, both in the upstream and downstream directions; therefore, 
no reflection of the reflected wave or of the transmitted wave by a 
terminal boundary is considered. 



determined exactly as before by considering the simultaneous solution 
of two boundary conditions at A: the continuity of pressure and the 

continuity of flow. The first boundary condition requires that 



The desired reflection and transmission coefficients will be 




( 1 . 21 ) 



and the latter requires that 



U 



A 




+ v' . 



( 1 . 22 ) 



Substituting Eqs. (1. 18) and (1. 20) into (1. 21) yields 




(1. 23) 



and substituting Eqs. (1. 17) and (1. 19) into (1. 22) yields 



(1. 24) 



1 1 



F 



A 




v 



Substituting Eq. (1. 15) for v and designating the resulting ratio by 



a = • Z , (1. 25) 

c 

the simultaneous solutions of Eqs. (1. 23) and (1. 24) result in the 
following values for the reflection and transmission coefficients, 
respectively: 

"' 26 > 



and 




2 g 

2 a + 1 



(1. 27) 



The sum of the squares of the two coefficients is not equal to unity 
as in the previous example by Lamb (1916), except when the value of 
a (and hence Z) is exactly equal to zero. The requirement that 
the value of Z given by Eq. (1.15) be equal to zero corresponds to 
no surface elevation in the resonator (e. g. , the resonator discharges 
into a large reservoir) or to an infinitely high velocity of discharge. 

The hydraulic impedance may be written in a form analogous to the 
electrical impedance; i. e. , 

Z = r + j (l CT ) ; (1. 28) 

Y U 

where r is proportional to the energy dissipative forces; and i , y 
are proportional to the inertial and potential (e. g. , gravity) forces. 
Requiring Z to vanish is analogous to a tuned resonant circuit in which 
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di = — — ; (1.29) 

Y ^ 

where (J is the angular frequency of oscillation. The value of Z will 
vanish at resonance only in an inviscid fluid (or in an analogous con- 
servative system). Equations (1. 26) and (1. 27) may be used for 
evaluating the amount of energy dissipation and the corresponding 
amount of reflection of waves propagating in a real fluid only when a 
value for r is known or may be determined. 

The difficulties encountered in the application of the results 
of Valembois (1953) lie in the evaluation of the parameter r and in 
the criteria for dimensioning a tuned resonator. The values of the 
hydraulic resistance and reactances are very small and are usually 
determined empirically. Secondly, an important omission in the 
theory is a rigorous analytical procedure for dimensioning the geo- 
metry of the resonator. The only reference given to a dimension was 
the empirically determined value for the length of the lateral expansion 
which was equal to one-quarter of a wavelength. Since surface gravity 
waves have no natural frequency of oscillation to which the natural 
frequency of the resonator may be matched, the electrical analogy 
cannot be employed as a means for dimensioning the resonator. 

This thesis is addressed to the evaluation of a reflection and 
of a transmission coefficient which are independent of any empirical 
data required for measuring the Small dissipative forces and which 
relate any geometric changes in the resonator to its effectiveness for 
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reflecting waves in terms of the characteristic parameters of the 
incident wavelength. 



2. THE SOLUTION 



Formulation of the Boundary Value Problem 



The derivation of the desired coefficients will proceed in a 
manner similar to that of the examples presented in Chapter 1 by the 
formulation of a boundary- value problem. A right-handed cartesian 
coordinate system (z-axis positive upwards) is oriented at the orthogo- 
nal intersection of the centerlines of the channel and of the resonator 
at the mean-water level (Figure 2. 1). The linear second-order homo- 
geneous hyperbolic partial differential equation which approximately 
describes the instantaneous surface elevation of the fluid in the domain 
mutually occupied by the channel and by the resonator is the wave 
equation; i. e. , 



2 2 
3 r\ _ 2 3 r| 

2 ~ c 2 

9 t 3 x 3 x 

k k 



( 2 . 1 ) 



where k has a range of values from unity to three. 

propagates at a constant celerity of 

O' 



The wave form 



( 2 . 2 ) 



Neglecting the non-linear effects of the boundary conditions requires 

that the following assumptions be made (McLellan, 1965): 

(1) The amplitude of the surface disturbance is 

small compared both to the wavelength and to 
the depth of the fluid. 
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Figure 2.1. The wave systems considered in evaluating the 
effect of the resonator on long waves. 
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(2) The channel is of uniform depth. 

(3) The fluid is inviscid and ir rotational. 

(4) The fluid is incompressible and homogeneous. 

(5) Coriolis acceleration may be neglected. 

(6) Surface tension may be neglected. 

(7) The rigid boundaries are smooth and impermeable. 

(8) The atmospheric pressure along the water surface 
is constant and uniform. 

The solution to Eq.(2. 1) for one dimensional space may be 
expressed in functional notation by 

t] = r| (x, t) . (2. 3) 

The most general solution is given by D'Alembert's Theory and is of 
the form 



r| (x, t) = F (kx - (F t) + G (kx + (Ft) ; 



(2. 4) 



where F is an arbitrary periodic surface disturbance propagating in 
the positive x direction and G is an arbitrary periodic surface dis- 
turbance propagating in the negative x direction. One simple -harmonic 
function which satisfies Eq. (2. 4) is a cosine function; i. e. , 



tl (x , t) = A, cos ( k x. dk (Ft + €. ) ; 
k j k j k 



(2. 5) 



where k has a range from unity to five corresponding to the com- 
ponents shown in Figure 2. 1 and j has a range from unity to two. 

The numerical indices for k are now replaced by alphabetical 
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indices shown in Figure 2. 1 in the following manner; 

1 = Mi 

2 = Mr 

3 = Mt 

4 = Ri 

5 = Er; 

where the upper case letters M, R correspond to the main channel 

and to the resonator components, respectively; and the lower case 

letters i, r, t correspond to the incident, reflected, and transmitted 

wave components, respectively. The components are replaced 

by x and y for the coordinate system shown in Figure 2. 1. Equation 

(2. 5) is a set of five wave component equations containing ten unknown 

constants; viz. , the five values for the wave amplitudes A,,., A,, , 

Mi Mr 

A,, , 'A_., A„ and their associated phase angles , 

Mt Ri Rr k s Mi Mr Mt 

e^., € . The set of equations expressed in Eq. (2. 5) are the following; 



r\... = A cos ( k x - CT t + e ) 
Mi Mi Mi 



r),. = A cos ( k x + (J t + e ) 

1 Mr Mr Mr 



•n _ , = A,, cos ( k x - dt + q, ) 
1 Mt Mt ' Mt 7 



11 Ri = A Ri cos < ky ' at + e Ri> 



•n = A^ cos ( k y + O’ t + e_ ). 
'Rr Rr ' y Rr ' 



( 2 . 6 : 1 ) 
( 2 . 6 : 2 ) 
(2. 6: 3) 
(2. 6:4) 
(2. 6: 5) 



In a well- stated boundard- value problem, the unknowns contained 
in Eqs. (2. 6; i) may be evaluated by a given set of boundary and initial 
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conditions (Tikhonov and Samarskii, 1963). In the stating of the con- 
ditions required by the above to comply with the criteria for a well- 
stated boundary-value problem, the following assumptions are made: 

(9) The fluid medium of propagation is at rest 
(i. e. , all velocities are absolute). 

(10) The terminal boundaries of the centerline of 
the main channel are sufficiently removed 
from the proximity of the resonator and, there- 
fore, are not required boundary conditions. 

(11) Steady state conditions exist. 

(12) All geometric dimensions considered are 
small compared to wavelengths (i. e. , long 
waves). 

(13) The cross sectional area of the channel and 
of the resonator consists of a horizontal 
bottom and parallel, vertical sides. 

(14) Shallow water conditions exist. 

Assumption (11) negates the requirement that the initial conditions 
be stated. Assumption (14) permits the use of the shallow water 



approximation for the wave celerity; i. e. , 



c = g h. (2. 7) 

The relative velocities shown in Figure 2.1 are related to Eqs. (2. 6 : i ) 
in the following manner; 



.A c 

77 Mi M ,, _ 

U = u * cos (k x - CT t; + e 

Mi h, . Mi 



U. 



Mr 



M 

^Mr C M 



M 



' cos ( k x + (Tt + € , ) 

Mr 



( 2 . 8 : 1 ) 
( 2 . 8 : 2 ) 
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U 

U 

U 



Mt 

Ri 

Rr 



A Mt C M 

h M 

A Ri C R 

h R 

A Rr C R 



R 



• cos (k x - dt + e Mt ) 

cos (k y - CTt + e ) 
i\l 



• cos ( k y + CTt + ). 

Rr 



( 2 . 8 : 3 ) 

(2. 8: 4) 

(2. 8: 5) 



The boundary conditions will be prescribed along the free surface in 
terms of Eqs. (2. 6: i ) and along the vertical boundaries in terms of 
Eqs. (2. 8: i ). 

The free surface boundary condition prescribes the instantaneous 
continuity of pressure at the origin for all time: 



11 Mi + 71 Mr " 11 Mt = 71 Ri + 11 Rr' 

The vertical boundary conditions require two statements, 
that the continuity of flow into tne domain is expressed by 



(2. 9) 

The first is 



- U w 1 b h + [U + U ] b h = 0; (2.10) 

Mi Mr Mt J MM L Ri Rr J R R 



and the second is that the no flow condition across the rigid terminal 
boundary of the resonator is given by 



[U . + U ] b_ h = 0. (2. 11) 

L Ri Rr J R R v ' 

The above three boundary conditions are evaluated at the following 

coordinates: 



Eq. (2. 9) at the origin (i. e. , x = y = 0); 

q _ b 

Eq. (2.10) at x = ±— , y = - -|i- ; 

and Eq. (2.11) at y = ~i , x = 0. 

R 
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Using the Boundary Conditions to Obtain a Solution 



The main channel incident wave may be arbitrarily selected as 
the reference wave profile by assumption (11) without any loss in 
generality; therefore, its associated phase angle is defined as 

‘Mi = °' < 2 - 12 > 

Eight independent boundary equations will result by the substitution of 
Eqs. (2. 6: i) and (2. 8: i ) into the three boundary conditions prescribed 
in Eqs. (2. 9), (2. 10) and (2. 11) when the expressions obtained after 
substitution are expanded by means of the trignometric identities for 
the sine and cosine of the sums and differences of angles. The time 
dependence of the boundary conditions may be eliminated by equating 
the sine and cosine terms which result from the above expansion. The 
ten equations which result from the preceding algebra are the following: 

(2. 13: 1) 
(2. 13: 2) 



A + A >, COS = COS e w._ 

Mi Mr Mr Mt Mt 



A, sine,, = A, , sine,, 
Mr Mr Mt Mt 



A,,. + A,, cos e,, = A„. cos e„. + A„ cos e„ 

Mi Mr Mr Ri Ri Rr Rr 



A,, sine,, = + A,, . sine^.-A,, sin e„ 
Mr Mr Ri Ri Rr Rr 



A,, cos e,, = A„. cos e_ . + A„ cos e„ 

Mt Mt Ri Ri Rr Rr 



A,, sin e,, = A^ . sin f . - A„ sin e„ 

Mt Mt Ri Ri Rr Rr 



(2. 13: 3) 
(2. 13: 4) 
(2. 13: 5) 
(2. 13: 6) 



A cos (k I - €_.) = A cos ( kl - e_ ) 

Ri ' R Ri' Rr ' R Rr 7 



(2. 13! 7) 
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-A,,, sin (k£_ -€_.) = A_, sin (k£ - €_ ) 

Ri ' R Ri 'Rr R Rr 



(2. 13: 8) 



b b b 

c.AJ A „;Cos(k-f-)-A^_cos(k— -^ Mr )-A Mt cos(k— + * W J] (2. 13:9) 



'M M L Mi 



Mr 



2 Mt' 



b b 

= c b [ -A . cos (k— €_,.) + A cos(k— e )] 

R R L Ri v 2 Ri' Rr ' 2 Rr J 



b b b 

-W + V Si ^T> + \!r Si " (k f-'Mr> +A Mt Sm(k f + 'Mt ' 1 < 2 ' 1 3! 10 » 



b b 

= c b_ [A^. sin(k - e^.) + A sin(k ^ -e )] . 

R R L Ri v 2 Rr Rr '2 Rr' J 



By means of Eqs. (2. 13: 7) and (2.1 3: 8), it may be shown that 



A 



Ri 




(2. 14) 



and 4 that 



+ 2k, R = 



Ri 



"Rr 



(2. 15) 



By squaring and adding Eqs. (2. 13: 5) and(2. 13: 6), Eq. (2. 14) may 
be rewritten in the form 



A 



Ri 



T2 



A 

Mt 

c o s ( ki R )[ ■ 



(2. 16) 



Equation (2. 7) may be employed to equate the two celerities; i. e. , 

(2. 17) 



c = c 
R M 



Following the substitution of the known relationships and applying some 
additional trignometric identities, the reflection and transmission 
coefficients may be obtained by squaring and adding Eqs. (2. 13: 9) and 
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(2. 13:10). The reflection coefficient may be expressed by 



b b 
R , M 



Sin (k ~2p + )[sin(k ^ ) - |tan(k^ R )| cos(k ™ )] 

M 



M 



K = 
R 



,l+(j^)sin(ky)[sin(k-y)- | tan(k£ R ) | cos^^)] + (^) 2 [ sin ( k_ ^ )' | ^(k^)) cos(ky-)] Z 
M M 



(2.18) 



and the transmission coefficient by 



k t = 



R 

cos (k 2 ) 



l + (r^)sin(k^)[sin(k^)- |tan(ktf R ) |cos(k ^)-] + ) 2 [sin(k“Y)-|tan(kjf R ) |cos(k ^)] 

- M M 



(2.19) 



The sum of the squares of Eqs. (2. 18) and (2. 19) is equal to 
unity and verifies that the energy in the system is conserved. 



Evaluating the Geometric Parameters 



The criteria for dimensioning the resonator are the maximizing 

and minimizing of Eqs. (2. 18) and (2. 19), respectively. This minimizing 

and maximizing process requires the evaluation of the arguments of the 

trignometric functions involving 0 , b . and b _ in terms of the tuned 

R R M 

wave length, and wave period, T Q , of the resonator. 

The dimension which has the most critical influence on the 
value of the coefficients is the transverse length of the resonator, j? , 
which appears in the argument for the tangent function. The tangent 
function becomes discontinuous when its argument approaches a value 



IV M 
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of (2n - 1 ) radians, where n is any positive mteger. The argu- 
ment reaches this value when 

L.. 

= (Zn " 1} ~T ; (2 * 20) 

where n must be equal to unity or two by assumption (12). Eq. (2. 20) 

agrees with the results obtained by Valembois (1953) when the value 

of one-half of the width of the main channel, b . is either a very small 

M 

per cent of or equal to one-half L^. The argument containing 

'i' 'r 

t may be expressed in terms of the angular frequency, O', by means 
ofThe wave celerity equation; i. e. , 



ki 



R 



fv 



Decreasing the dimension H in Eq. (2. 21) to a value less than 

R 



( 2 . 21 ) 

L * 



increases the tuned resonant frequency, (J , by a corresponding 

° L * 

decreases 



factor; while increasing 9 to a value greater than — ■ 

R 4 



the tuned resonant frequency, CT , also by a corresponding factor; i. e. , 

L.. ° 

if i R < 



— — , then 0* > O' 
4 o 



L^ 

and if i ^ > — - — , then 0*< O' ; 

where O’ is the tuned resonant frequency for maximum reflection. It 
may be seen from Eq. (2. 16) that the resonator wave amplitude would 
become infinite at resonance for any finite value of the transmitted 
wave height, A . Since an infinite wave amplitude does not 
seem to be physically valid, the transmitted wave height, 
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must be identically zero. 

The influence of the width of the resonator, b , on the value 

R 

of the coefficients is most easily evaluated by means of Eq. (2. 19). 
Minimizing the transmission coefficient requires that the value of the 
cosine function in the numerator approach zero. The cosine function 
approaches zero when its argument approaches (2n -1)*—- , where n 
is any positive integer. Therefore, the value of the width of the 
resonator must be 



b R = Un-1>— ! 



( 2 . 22 ) 



where n must be unity by assumption (12). If assumption (12) is 

valid, the value of b determined by Eq. (2. 22) is a maximum; 

R 

hence, b may only decrease from this maximum value. The argu- 

JL\ 

ment containing b may be expressed in terms of the angular frequency 

R 



as before by 



R 



O' b R 



c 2 



(2. 23) 



Decreasing b from its maximum value given by Eq. (2. 22) increases 
R 

the tuned resonant frequency CT . 

o 

Finally, ''a critical value for the width of the main channel, b , 

M 

is not easily evaluated by means of the equations for the coefficients; 
however, because the preceding derivation was based on the assumption 
that all -dimensions were small compared to a wavelength, it may be 
determined that the width of the main channel must be a small per cent 
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of the tuned wavelength, L , . 

The evaluation of the geometry of the resonator system has 
been discussed by considering the effect of a geometrical change on 
the ability of the resonator to reflect an incident wave of a single fre- 
quency. Assumption (11) does not allow the application of the coeffi- 
cients to a constantly varying input inducing a transitory response. 
However, the experimental results obtained by Valembois (1953) 
indicate that appreciable reflection does occur over a range of fre- 
quencies. In fact, the response curves obtained by Valembois are 
similar to the response curves for electrical filter circuits (Olson, 
1958). The analogy of these two oscillatory systems introduces the 
possibility of determining an effective frequency band-width for which 
reflection may be significant. 



3. PRESENTATION OF WAVE DATA 



The derivation of the reflection and transmission coefficients 
in Chapter 2 and their subsequent discussions dealt with an incident 
wave of a single frequency and wavelength. An analysis of an incident 
wave record reveals that waves in nature are composed of an infinite 
and continuous number of frequencies resulting from a complex 
phenomenon. Since the resonant structure is static, some means 
of determining a single frequency which is a characteristic measure 
of the complicated group phenomenon contained in the waves must be 
made as well as some means for evaluating the effect of the resonator 
when subjected to an input composed of an infinite series of incident 
frequencies. Fortunately, the solution to the first problem provides 
a means for resolving the second. 

The problems encountered in recording and in analyzing surface 
waves are currently under research by the discipline of Oceanography. 

A complete or even adequate treatment of these complex problems is 
beyond the scope of this thesis. The design and construction of the 
resonant breakwater serves to illustrate the interdependance required 
between Civil Engineering and Oceanography in oceanographic con- 
struction. The following discussion is intended only to provide some 
background and continuity for the use of the oceanographic data required 
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to effect a complete solution to the design problem posed. The refer- 
ences found in the Bibliography will provide a means for initiating a 
more detailed study concerning spectral analysis of ocean waves. 

Recording Wave Data 

The answer to the rhetorical question of how to extract a 
characteristic frequency from those occurring in a natural phenomenon 
is affected by the methods available for recording the actual physical 
profile of the surface disturbance. The various recording devices 
currently in use all have the common failure of a lack of a standard 
measurement and of distortions inherent in the recording apparatus. 

The effect of the recording distortions becomes more pronounced if 
small non-linear components must be recorded accurately; but the 
presence of these inherent distortions will influence the shape of any 
wave profile recorded and will, therefore, be incorporated in the data 
for a nonlinear analysis as well. Because the recording process is 
quite complicated, the author can only refer to oceanographic literature 
and follow Kinsman (1965) by assuming that for the case in question 
an adequate and representative wave record is available to the 
designer from the oceanographer. 
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The Wave Spectrum 

The assumed wave record obtained from the oceanographer 
may also provide an answer to the again rhetorical question of how to 
extrapolate a characteristic design frequency from the natural wave 
phenomenon. The continuous distribution of the frequencies in a wave 
train suggest that a statistical technique may be employed to extract 
the characteristic frequency. The Sverdrup and Munk Theory, for 
example, describes the sea surface by a single sinusoidal component 
called the "significant wave" (Kinsman, 1965). The next level of 
complexity involves a Fourier analysis. 

In a Fourier or harmonic analysis, the wave record is mechani 
cally transformed into a spectrogram (McLellan, 1965). If the charac 
teristic solutions to Eq. (2. 1) are orthogonal with respect to the 
weighting function (which arises from the Sturm - Ldouville problem 
generated by the separation of variables) and are piecewise differenti- 
able in the recording interval T (Hildebrand, 1962), the solutions 
may be represented by a product series of the form 

oo 

T|(t) = / t,. (t). (3. 1) 

L-j n 

n = 1 

The solution determined in Chapter 2 was a cosine function. Since 
this function is symmetrical with respect to the z-axis, it may be 
expressed in a Fourier series for an even function; i. e. , 
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oo 

V- H 

n (t) = ) COS ( <J t - € ) . (3. Z) 

n L-j c n n 

n = 1 

The coefficients H are determined by multiplying each side of Eq. 

(3. 2) by the weighting function and by the k u solution, , given 
in Eq. 3. 1, and then integrating over the truncated recording interval, 

T. 

00 

H = ~ \ T| (t) • cos (O' t - e ). (3.3) 

n T v q n n 

The value given by Eq. (3. 3) is equal to twice the average 
value of the product r\ (t) * cos ( CT t - € ) over the recording interval 

T. The wave energy spectrum may also be derived from Eq. (3. 2) 
by squaring both sides of the equation and integrating over the recording 
interval T. The harmonic analysis of the wave record does not pro- 
vide any information regarding the shape of the surface disturbance 
when 

| t \ > T (Kinsman, 1965) . (3.4) 

The next logical progression from this discrete spectrum is a repre- 
sentation of the wave record by a continuous spectrum. 

The Bretschneider Continuous Spectrum 

The purpose of this discussion on the methods of presentation 
of wave data is to enable the designer either to .request wave data in 
a particular form or to transform a wave record available from the 
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oceanographer. One such method of wave data presentation currently 
in use is the Bretschneider Power Spectrum. This spectrum requires 
that the average height of the highest one-third waves, H^^y as 
determined from a Rayleigh probability distribution of wave heights, 
be known (Dean, 1966). From the value of the significant wave height, 
iy a parameter which is equal to twice the total energy contained in 
the waves is defined by 

H 1/ , 3 = 2. 83 \T E f . (3. 5) 



Bretschneider determined that the spectral distribution of the wave 

energy, S 2» as a function of frequency CT could be represented by 
H 



v 1 "- 






CT 



0 “ 



(3.6) 



The parameter 0" is the frequency at which the peak energy occurs 



and is defined in terms of the significant wave period, T^^» ky 



CT 

o 



2 t r 



1. 17 




(3.7) 



Integrating Eq. (3. 6) from zero to infinity yields twice the energy 
contained by the waves E^. 



A Functional Design 



The design of the resonant breakwater considered in this thesis 
is an example of a functional design discussed by Wiegel (1964); 
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however, only part of the protection function (which is but one aspect 
of the over-all function of the breakwater) is developed. The functional 
criteria also affect the method required for presenting the wave data. 
The purpose for which the harbor is being constructed dictates both the 
functional requirements of the resonator and the required method of 
wave data presentation. 

As discussed previously, the wave height and wave period 
parameters normally available from the oceanographer are those of 
the significant wave, and T^^y The design criteria will 

normally require protection from the effects of the highest waves 
present at the structure. The significant wave height, has 

been determined to be an accurate index to structural damage by 
waves (Wiegel, 1964). 

Harbor oscillations or seiching will also normally be found 
in the design criteria for a harbor. Seiching is a complicated response 
of the harbor to a periodic force. This response is a function both 
of the geometry of the harbor and of the period of the force. Wave 
data for designing protection from seiching, therefore, must be a 
function of frequency. 

One solution to the problem posed at the beginning of this 
chapter of how to determine a single design frequency from the con- 
tinuous spectrum of frequencies occurring in the natural phenomenon 
appears to lie with the statistical significant wave period, T^^y 
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The Bretschneider Power Spectrum derived from the significant wave 
height, will then provide a means of evaluating the effectiveness 

of the resonator to an input composed of an infinite and continuous 
spectrum of frequencies. By means of this power spectrum, the effect 
of a resonator on a wave system occurring in nature may now be 



evaluated. 



4. AN EXAMPLE OF A FUNCTIONAL DESIGN 



The Bretschneider Power Spectrum offers a solution to the 
design problems posed at the beginning of Chapter 3. Evaluating the 
effectiveness of the resonator over a spectrum of frequencies by com- 
bining this power spectrum with the reflection and transmission coef- 
ficients represents a linear analysis. The analysis is linear since 
both the power spectrum and the expressions for the coefficients were 
derived from linear solutions. 



A Linear Analysis 



In a linear process, the output, 0, is a linear function of the 
input, I ; i. e. , 



0 = R (0") • I . (4. 1) 

R(CT) is the linear transfer function and is a function of frequency, (T. 
The essential features of a linear process are the following (Dean, 1966) 

1. Changing the amplitude of the input function, 

I, changes the amplitude of the output function, 

0, proportionally. 

Z. The constant of proportionality is the linear 

transfer function and is a function of frequency, 

(7, only (i. e. , the transfer function is inde- 
pendent of the amplitude). 

3. Input functions may be superimposed. 

The effectiveness of the resonator at a given frequency, therefore, 
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may be determined from the product of the square of Eq. (2. 18) or 

(2. 19) and the value of the ordinate S -> from Eq. (3. 6) for the 

n 

Mi 

same corresponding frequency; i. e. , 

s 2 (CT ) = (<T ) • S 2 (O';) ( 4 - 2 ) 

H 1 K 1 H 1 

Mr Mi 

or 

s 2 (°\> = K 2 - ( O’ ) • S 2 (O'. )• (4.3) 

1 T 1 H 1 

Mt Mi 

Wave Data and a Linear Analysis 

As an illustrative example of a linear analysis, it is assumed 
that a previously undeveloped coastal area has been selected for the 
construction of a new harbor facility. The harbor has been designed 
to accomodate large commercial vessels (Hennes and Ekse, 1955). The 
proposed entrance to the harbor between the breakwaters is to have 
a depth of fifty (50) feet and a width at the mean-water level of two- 
hundred (200) feet (Figure 4. 1). The effectiveness of the construction 
of a single resonator on the local wave conditions is to be measured. 

The preliminary engineering report includes the collection of 
extensive oceanographic data on the proposed location. From this 
report, the significant wave height and period were determined to be 
the following: 




o 
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Figure 4.1. Proposed harbor entronce design. 
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H^ 3 = 1 2 feet 

= 1 2. 7 sec 

The total energy contained in the waves is determined from Eq. (3. 5) 
to be 

= 18 sq. ft. 

The significant wave period, T^^, determines the value of CT^ in 

Eq. (3.7) from which the spectral distribution of the incident wave 

energy, S ^ , may be determined by means of Eq. (3. 6) (Figures 

H Mi 

4. 2 and 4. 3). 

By satisfying the requirement of Eq. (1. 6), the shallow water 
approximation for the wave celerity is determined to be forty (40) feet 
per second through the breakwater entrance. By means of Eq. (2. 2) 
and the value of the wave celerity, the significant wave length, L» f , is 
found to be six-hundred (600) feet. Utilizing the design parameters 
from Eq. s (2. 20) and (2. 22), the dimensions of the resonator are as 
follows: 

f R = 150 feet 
b R = 300 feet 

The values of the reflection and transmission coefficients may now 
be determined as a function of frequency from Eqs. (2. 18) and (2. 19), 
respectively. The frequency dependence of the square of these 
coefficients is also shown in Figures 4. 2 and 4. 3. 
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T (sec) 





MAIN CHANNEL INCIDENT WAVE ENERGY SPECTRUM 




MAIN CHANNEL REFLECTED WAVE ENERGY SPECTRUM 
Figure 4.2. Determination of the effective frequency band-width of the resonator. 
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T (sec) 




SQUARE OF THE TRANSMISSION COEFFICIENT 




MAIN CHANNEL INCIDENT WAVE ENERGY SPECTRUM 




MAIN CHANNEL TRANSMITTED WAVE ENERGY SPECTRUM 
Figure 4.3. Determination of the effective frequency band-width of the resonator. 
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The reflective and transmissive effectiveness of the resonator 
at a given frequency is determined from the product of the ordinate of 
the incident wave spectrum and the ordinate of the transfer coefficient. 
The resulting spectral distributions of these two linear analyses are 
shown on Figures 4. 2 and 4. 3. 

Discussion 

The resonator appears to have a frequency band- width ranging 

from 13 to 17 seconds in which reflection is approximately 100%. The 

effectiveness of the resonator is determined from the amount of energy 

which it reflects and transmits. This amount of energy is measured 

from the area under the power spectrum curves for the reflected and 

transmitted waves, S 2 and S 2 > shown in Figures 4. 2 and 

^ Mr H 

4.3, respectively. The narrow spike occurring between the periods 
of 13 and 14 seconds prevents the resonator from being greater than 
85% effective in reflecting the incident wave energy. The proposed 
design for the given wave conditions is approximately 75% effective 
for reflecting the energy contained in the incident wave spectrum 



between the periods from 10 to 17 seconds. 
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